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Search for new physics in tt + ^t -^ bbqqqq + ^t final state in pp colhsions at 

Vs = 1.96 TeV 
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We present a search for a new particle T' decaying to a top quark via T' — > i + X, where X goes 
undetected. We use a data sample corresponding to 5.7 fb~^ of integrated luminosity of pp collisions 
with y/s — 1.96 TeV, collected at Fermilab Tevatron by the CDF II detector. Our search for pair 
production of T' is focused on the hadronic decay channel, pp — >■ T'T' — >■ tt-\- XX — >■ bqqbqq -\- XX. 
We interpret our results in terms of a model where T' is an exotic fourth generation quark and X is 
a dark matter particle. The data are consistent with standard model expectations. We set a limit 
on the generic production of T'T' -^ ti -\- XX, excluding the fourth generation exotic quarks T' at 
95% confidence level up to my/ = 400 GeV/c^ for mx < 70 GeV/c^. 
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There are many hints, from astronomical observations 
and cosmological theories, for the existence of dark mat- 
ter (DM) particles, which must be long-lived on cosmo- 
logical time scale [l|. The long lifetime of DM can be 
explained by the conservation of a charge of a new sym- 
metry. Direct-detection experiments based on ultra-low 
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noise devices have recently obtained interesting results. 
The DAMA/LIBRA Collaboration (3], searching for an- 
nual modulation in the interaction rate due to the earth 
motion through DM galactic halo, has claimed a ~ 9(7 
observation of DM. The CoGeNT Collaboration has also 
reported evidence of DM [3|. If confirmed, these results 
would imply, unlike astronomical observations, that DM 
interactions with standard model (SM) particles are not 
only gravitational. DM interactions with SM particles 
could be allowed by weak interactions, or by connector 
particles carrying both dark and SM charges, so that 
they could be produced in colliders. The second hypoth- 
esis is favored in the case that DM particles have mass 
of a few GeV/c^, as DAM A and CoGeNT results seems 
to indicate. In a recent model [4] the role of a connec- 
tor particle is played by an exotic fourth generation T' 
quark, which is supposed to decay to a top quark and 
dark matter, T' ^ t + X . The pair production of such 
exotic quarks and their subsequent decay has a collider 
signature consisting of top quark pairs (it) and missing 
transverse energy (^t) @ due to the invisible dark mat- 
ter particles. These types of signals are of great interest 
as they appear also in other models containing DM can- 
didates, such as scalar top quarks production and their 
decay to top quarks and neutralinos p| or top quarks 
and gravitinos [7], and in many other new physics sce- 
narios such as little Higgs [8| and models where baryon 
and lepton numbers represent local gauge symmetries [9j . 

A first search for the T'T' -^ ti + XX process has 
been performed in the semileptonic channel: ti+XX — >■ 
bWbW + XX -> bliybqq + XX ^. This letter reports 
the first search for such a process in the all-hadronic tt 
decay channel, characterized by a larger branching ra- 
tio and a lower physics background rate. Events were 
recorded by CDF II ^T] a general purpose detector de- 
signed to study collisions at the Fermilab Tevatron pp 
collider at ^/s = 1.96 TeV. The tracking system consists 
of a cylindrical open-cell drift chamber and silicon mi- 
crostrip detectors in a 1.4 T magnetic field parallel to the 
beam axis. Electromagnetic and hadronic calorimeters 
surrounding the tracking system measure particle ener- 
gies and drift chambers located outside the calorimeters 
detect muons. Jets are reconstructed in the calorime- 
ter using the jetclu [Ij] algorithm with a clustering 
radius of 0.4 in azimuth-pseudorapidity space [1J|- The 
detector response for all simulated samples is modeled 
by a detailed CDF detector simulation. Production of 
T' pairs and their subsequent decays to top quark pairs 
and two dark matter particles would appear as events 
with missing transverse energy from the two dark mat- 
ter particles, and six jets from the two b quarks and the 
hadronic decays of the two W bosons. We model the 
production and decay of T' pairs with the MADGRAPH 
Monte Carlo (MC) generator [IJ], and normalize to the 
next-to-next-to-leading order (NNLO) cross section cal- 
culation [15|. Additional radiation, hadronization and 
showering are described by PYTHIA |16j . 

We use a data sample corresponding to an integrated 



luminosity of pp collisions of 5.7 fb~ , collected by re- 
quiring ^T > 50 GeV and two or more jets with trans- 
verse energy Et > 30 GeV and \rj\ < 2.4 . We then 
require 5 < Nj^ts < 10, where Njets is the number of 
jets, and where all jets satisfy the requirement \ri\ < 2.4. 
We also require the transverse energy E^ of the sublead- 
ing jets, Ji, to be greater than 20 GeV for (i = 3,4,5) 
and 15 GeV for (i > 5). We veto events with at least 
one isolated electron or muon to suppress events with 
semileptonic ti decay. We refer to this sample as the 
preselection sample. At this stage of the event selec- 
tion, multijet QCD background where ^t arises from jet 
energy mismeasurement accounts for more than 95% of 
the expected backgrounds. The second dominant back- 
ground is tt production. We model this process using 
PYTHIA with mt = 172.5 GeV/c^ [17], normalized to the 
next-to-leading (NLO) order cross section [1^. Associ- 
ated production of W/Z boson and jets is also a signifi- 
cant background source. Samples of simulated W/Z+jets 
events with light- and heavy-flavor jets are generated us- 
ing the ALPGEN [1^ MC generator, interfaced with the 
parton shower model of PYTHIA. A matching scheme is 
applied to avoid double-counting of partonic event con- 
figurations [20| . The W/Z +]ets samples are normalized 
to the measured W and Z cross section [2l[. Diboson 
and single top production are modeled using respectively 
PYTHIA and MADGRAPH, and normalized to NLO cross 
sections |2l| - |24| . Because of the large production rate 
for QCD multijet events at a hadron collider and the 
statistics needed in order to describe this process ade- 
quately in an analysis looking for a very small signal, the 
Monte Carlo simulation of QCD multijet events is pro- 
hibitive. More importantly, the systematic uncertainties 
associated with the Monte Carlo simulation of QCD jet 
production are large. For these reasons, we estimate the 
QCD background solely from data. Similarly to Et, it 
is possible to define a missing transverse momentum pT 
using the spectrometer, as the negative vector sum of the 
charged particles momenta. Et and ^t are correlated in 
magnitude and direction in events with undetected par- 
ticles. In QCD multijet events Et originates from the 
mismeasurement of a jet energy in the calorimeter, while 
Pt depends on fiuctuations in the number of charged par- 
ticles in a jet, so they are usually aligned or anti-aligned 
in dijet-like events like energetic QCD multijet events, as 
is shown in Fig. [TJ QCD multijet events in which Et 
and "^T are aligned or anti-aligned have the same kine- 
matic characteristics, as we have verified studying QCD 
multijet samples with 2 and 3 jets [2^. We reject events 
with A<I){Et,^t) > 7''/2, and use them to model QCD 
multijet events in the signal region A(f>{ET,^T) < 7r/2. 
To further suppress the QCD multijet background, we 
require the azimuthal distance between the directions of 
Et and subleading jets, A(J){Et, Ji), to be greater than 
0.4 for i = 1,2,3 and 0.2 for i = 4, 5. We also require 
^T > 20 GeV and ^t sig > dGeV^/"^, where ^t sig is 
defined as the ^t divided by the square root of the total 
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FIG. 1: Distribution of A0(.^t,^t) for the preselection 

data, and two scenarios with different values of rriT' and 

TTT-x- AH histograms are normalized to unit area. 



energy collected in the calorimeter. Finally, we require 
J2jets ^T > 220 GeV to remove soft QCD events. Ah 
these cuts have been chosen to optimize the 5/^/(5 + B) 
figure of merit, where S and B are respectively the ex- 
pected numbers of signal and background events. TableU 
shows the expected number of events in the signal region 
for SM backgrounds and for several signal hypothesis. 



TABLE I: Number of expected signal events for five 

benchmark scenarios compared to data and expected 

SM backgrounds. 






ttXX(hadronic) [GeV/c^ 



Events 



niT' 


, mx = 260, 80 


88.5 ±11.9 


rriT' 


mx = 330, 100 


66.4 ±8.9 


rriT' 


mx = 360, 100 


39.7 ±5.3 


rriT 


>,mx = 380,1 


27.3 ± 3.7 


rriT 


, mx = 400, 1 


17.5 ±2.3 




QCD 


745.4 ± 124.3 




tt 


498.2 ± 66.8 




W+jets 


119.7 ±48.4 




Z+jets 


39.4 ±15.9 




Diboson 


17.9 ± 2.2 




Single top 


5.3 ±0.8 



Total Background 



1423 ± 150 



Data 



1507 



Inverting one of the event selection cuts, keeping oth- 
ers unchanged, allows us to define a signal-depleted con- 
trol region. We use ^t sig < 3GeF^/^, Njets = 4 and 
•^T < 20 GeV control regions to validate the overall back- 
ground modeling. The normalization factor of the QCD 
background is given by the average ratio of QCD events 
that pass the A(J){Et,^t) < 7r/2 cut to QCD events that 
fail the cut in these three control regions. Figure [H shows 
good agreement of background modeling with data in 
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FIG. 2: Top plot shows the ^t sig distribution in 

events with four jets and large ^t- Bottom plot shows 

the Njets distribution in event with 5 < Nj^ts < 10 and 

^T sig < 3 VGeV- 



these regions. We consider several sources of systematic 
uncertainties. The dominant components are the uncer- 
tainties on the QCD normalization factor, the jet en- 
ergy scale (JES) ^] and the theoretical cross sections. 
We also take in account the differences of ti predicted 
rates using different hadron fragmentation models in the 
HERWIG [27] Monte Carlo, and varying initial/final state 
radiation and color reconnection effects [2g] . The vari- 
ation of the JES was found to change significantly the 
^T sig distribution in addition to its normalization, and 
its variation is thus taken into account. Figure [3] shows 
the ^T sig distribution for expected signal events and 
SM backgrounds. The signal is expected to contribute 
significantly in the high tail of the ^t sig distribution. 
There is no evidence for the presence of T' — )■ t + X 
events in the data. We calculate 95% C.L. upper limits 
on the T' ^ t + X cross section, by performing a binned 
maximum-likelihood fit on the ^t sig distribution. The 
limits are calculated using a Bayesian likelihood method 
with a flat prior for the signal cross-section, integrat- 
ing over Gaussian priors for the systematic uncertain- 
ties. The results are shown in Table |lTl We convert 
the observed upper limits on the pair-production cross 



sections to an exclusion curve in mass parameters space 
{mT',mx)- As shown in Fig. SI a significant enhance- 
ment in sensitivity is obtained when comparing to the 
previous analysis in semi-leptonic channel. 



exclusion range up to m^' = 400 GeV/c^, for mx < 70 
GeV/c^ . Finally, this study shows that the bbqqqq + ^t 
final state is the most sensitive to the generic production 
of top quarks plus dark matter candidates, and thus the 



TT' ^ tt + XX 

(m^,, m^) = (330, 40) GeV/c^ 

(m^,, m^) = (380, 1 ) GeV/c^ 
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FIG. 3: ^T sig distributions for the standard model 

backgrounds, the observed data, and for two scenarios 

with different values of uit' and mx ■ 




mr=m,op+mx 
expblvbqq+XX4.8fb' 
obsblvbqq+XX4.8fb' 
exp bqqbqq+XX 5.7fb I 
obs bqqbqq+XX 5.7fb^ 



200 



250 



300 



350 



400 



mp (GeV/c') 



TABLE II: Expected 95% C.L. upper hmits on T'T' 

production, where the uncertainty range covers 68% of 

the pseudoexperiments, and observed exclusion limits 

for representative signal points. 



{mT',mX) GeV/c' 


(^exp.95%C.L.excl. {pb) (^obs.95%C.L .excl.{pb) 


(200,40) 


2.02 ± 0.65 1.90 


(220,40) 


2.14 ± 0.75 3.00 


(260,1) 


0.23 ± 0.08 0.18 


(280,1) 


0.15 ± 0.05 0.12 


(280,40) 


0.18 ± 0.07 0.15 


(300,1) 


0.09 ± 0.03 0.09 


(300,80) 


0.20 ± 0.06 0.16 


(300,100) 


0.29 ± 0.09 0.38 


(330,1) 


0.05 ± 0.02 0.03 


(330,100) 


0.13 ± 0.04 0.18 


(360,1) 


0.03 ± 0.01 0.02 


(360,100) 


0.06 ± 0.02 0.04 


(380,100) 


0.06 ± 0.02 0.05 


(400,1) 


0.023 ± 0.008 0.016 



In conclusion, we performed the first search for new 
physics in the tt + ^r -^ bbqqqq + 6t final state. Data is 
consistent with the background-only hypothesis, and we 
thus set 95% C.L. upper limit on the production cross 
section for fermionic T' pairs decaying to top quarks and 
dark matter candidates X, increasing the existing mass 



FIG. 4: Expected (exp) and observed (obs) 95% C.L. 
exclusion region in the {mT' ,mx) parameters space. 



most promising to probe the supersymmetric i — >■ t -I- xl 9 
scenarios at the LHC. 
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